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Abstract. 

High resolution observations of ^■^CO and ^^CO 1-0 in the Medusa 
(NGC 4194) minor merger show the tj§§ 1~0 intensity ratio (7^i_o) in- 
creasing from normal values (5-10) in the outer parts of the galaxy to 
high (> 20) values in the central, extended starburst region. Ratios > 20 
are otherwise typical of more luminous mergers. The Medusa ^™ ratio 
rivals that of ultraluminous galaxies (ULIRGs), despite the comparatively 
modest luminosity, indicating an exceptionally high star formation effi- 
ciency (SFE). We present models of the high pressure ISM in a ULIRG 
and the relatively low pressure ISM of the Medusa. We discuss how these 
models may explain large 7^i_o in both types of distributions. Since the 
HON emission is faint towards the Medusa, we suggest that the SFE is 
not primarily controlled by the mass fraction of dense (nii lO^cm^'^) gas, 
but is probably strongly dependent on dynamics. The bright HCN emis- 
sion towards ULIRGs is not necessarily evidence that the IR emission 
there is always powered by starbursts. 



1. Introduction 

Infrared luminous (Lir > 10^^ Lq) and ultraluminous galaxies (ULIRGs) (Lir > 
10^^ Lq) are believed to be major mergers of massive, gas-rich disk galaxies. 
They are spectacular systems, related to super-starbursts and massive gas flows 
to their central regions. Their molecular gas distribution is often found in ro- 
tating nuclear disks, toruses or bars only a few hundred pc in extent. These 
ULIRGs have been suggested as precursors of QSOs (e.g. Norman and Scoville 
1988). Still, galaxies collide that do not have the necessary properties to become 
ultraluminous objects. The importance of these intermediate, often minor i.e. 
unequal-mass mergers, to the general evolution of galaxies is not well under- 
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Figure 1. The left panel shows an overlay of the ^^CO (not primary 
beam corrected) contours (white, apart from clouds in the tail marked 
in black) on a grcyscalc, overexposed optical R-band image (Mazzarella 
Sz Boroson 1993). The dashed curve marks the edge of the OVRO 
primary beam. The right panel shows the CO velocity field. Contours 
range from 2400 to 2640 kms~^, with spacing 24 kms~^. The grayscale 
ranges from 2400 (light) to 2650 (dark) kms~^. The radio continuum 
peak is marked with a cross. Note that the position angle of the ^^CO 
distribution is perpendicular to the major axis of the velocity field 
because of the crossing of the central dust lane. 



stood. Both major and minor mergers undergo bursts of star formation, but 
the activity seems to take place on a larger linear scale in the minor mergers. 
Whether the starburst processes and triggering mechanisms remain the same for 
the compact and large scale bursts is unclear. It appears however, that some of 
the extended bursts can reach the same high SFEs as the ULIRGs. 

Probing the properties and distribution of the molecular ISM enables mod- 
elling of the triggering and evolution of the starburst. For this purpose, it is 
necessary to go beyond ^^CO and use other (fainter) molecular tracers such as 
the isotopomer ^^CO and high dipole moment molecules such as HCN, HNC 
and CN. These lines are powerful diagnostic tools, in particular in combination 
with other extinction-free tracers like IR and radio. 

The minor merger, and young shell galaxy, NGC 4194 ("the Medusa") has 
an extended region (2 kpc) of intense star formation (e.g. Prestwich et al. 1994) 
which is responsible for most of the FIR luminosity (-Ljr = 8.5 x 10^*^ Lq at 
L»=39 Mpc). The SFE, is high, 40 Lq Mq\ similar to the SFE of the 

ULIRG Arp 220. The global TZi-q is large, « 20, indicating a highly excited or 
disrupted ISM (Aalto &; Hiittemeister 2000, AH). Despite the high SFE, emission 
from the high density tracer molecule HCN remains undetected in NGC 4194. 
We present a preliminary analysis of the ^^CO and ^'^CO distribution obtained 
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at high resolution with OVRO. We offer an explanation for the apparent lack 
of dense gas in the Medusa and suggest "cartoon models" of a high pressure 
(ULIRG) and a comparatively low pressure (Medusa) molecular ISM. 



2. The isqq 10 intensity ratios 

Since both high kinetic temperatures and large turbulent line widths will de- 
crease the ^^CO 1-0 optical depth (tco), a map of the n^Q 1-0 intensity ratio 
(7^i_o) can be used to identify regions of extreme or unusual physical conditions 
in the molecular gas. Large ratios indicate low to moderate rco a-^d Aalto ct al. 
(1995) established some general diagnostics of the cloud conditions and environ- 
ment based on global values of TZi-q: Small ratios, 7?.i_o~ 6 are an indication of 
a normal Galactic disk population of clouds dominated by cool giant molecular 
clouds (GMCs); intermediate ratios 10;:7^i_o^l5 are associated with the inner 
regions of normal starburst galaxies; the extreme values TZi-o > 20 originate 
in turbulent, high pressure gas in the centers of luminous mergers. In the most 
luminous mergers, the gas surface density implied by the Galactic conversion 
factor from -"^^CO luminosity to M(H2) is well over 10^ Mq pc~^ — two orders of 
magnitude higher than in typical Milky Way GMCs. This led Aalto et al. (1995) 
to suggest that large values of TZi-q are related to the large gas surface densities 
in compact nuclear starbursts. Large surface densities require high pressures in 
hydrostatic equilibrium; as the ^^CO line ratios in these objects indicate 
low {n < 10^ cm~^) densities of the ^^CO emitting gas-component it must be 
supported by large turbulent line widths {P oc n{6V)'^). Thus, rco can be re- 
duced to moderate (r^ 1) values, resulting in large TZi-o- However, the global 
7^1-0 of the Medusa is high, but the gas surface density is more than an order of 
magnitude lower than that of the ULIRGs. Therefore, the pressure will be con- 
siderably lower and the ISM is likely much less turbulent (even though it is still 
significantly higher than the average pressure for the Galactic disk). Thus, for 
the Medusa, the high 7^i_o niay well be caused by elevated kinetic temperatures 
rather than large line widths. There is indeed a strong correlation between FIR 
colour temperature and large values of TZi-q: Galaxies with 7^i-o > 20 all have 
/(60/xm)//(100/xm) flux ratios >0.8 (e.g. Aalto et al. 1991, 1995) indicating 
high average dust temperatures. Thus, large values of TZi-o are also related to 
the effect of the starburst itself, heating dust and gas to high temperatures. 

2.1. i2co and ^^CO morphology of NGC 4194 

b igure 1 shows the ^^CQ emission overlayed on an R-band image of NGC 4194. 
The ^^CO distribution is surprisingly extended, 5 kpc, for an advanced merger 
even if the bulk (60 - 70 %) of the emission emerges from the 2 kpc starburst 
region. The morphology is complex, tracing large-scale dust lanes, one which 
curves along the north-eastern edge of the main body and continues into the 
tidal tail, and one which is crossing the central region along the minor axis. The 
gas mass (for a standard CO to H2 conversion factor) is 2 x 10^ Mq (AH). There 
is substantial variation in TZi-o across NGC 4194: From quiescent values of 7 
in the the eastern dust lane 5 kpc from the center, to high values, >20, in the 
starburst and central dust lane. This clearly demonstrates that there is a strong 
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Figure 2. Left image shows CO contours overlayed on HST WFPC2 
image, and the right shows the ^^CO 1-0 contours. The contours for 
both images have levels (20,40,60,80,90)% of the peak flux. The reso- 
lution is 4" for both CO and ^^CO. The dashed line outlines the dust 
lane and the cross the radio continuum peak. 



connection between 7^i-o ^-^d gas environment. In Figure 2 we show the CO 
and ^^CO contours overlayed on an archival HST WFPC2 image of the inner 2 
kpc of NGC 4194. The ^^CO and ^^CO morphologies differ substantially showing 
significant variation in 7^i_o also within the central region. ^^CO peaks west of 
^^CO, associated with a peculiar dust feature. Elevated values of 7^i-o are found 
in two regions: a) througout the extended starburst region and b) in the curved 
part of the central dust lane. The global TZi-o is dominated by the starburst 
region where 7^i-o is likely elevated because of large kinetic temperatures (see 
next section). In region b) gas streaming may cause local effects of large line 
widths (see also Hiittemeister and Aalto, this volume). 



Table 1. Line Ratios in high 7^i-o starbursts 

Line High pressure (compact) Low pressure (extended) 
CO 2-1/1-0 subthermal (<0.6) thermahzed (« 1) 

7^l_o ^20 >20 

CO/HCN 8 >25 



2.2. The molecular ISM of an extended and compact starburst 

In the table we list typical line ratios for compact high pressure and extended 
low pressure starbursts. Both types may have large values of TZi-o, but their 
average ISM properties differ on several accounts. In Figure 3 we show cartoons 
of the two ISM types, "raisin roll" and "fried eggs" . For the former (ULIRG) 
scenario, the ^^CO emission is emerging from low density n^lO^ cm""^ diffuse 
gas of large linewidths and filling factor, while the HCN 1-0 emission is coming 
from dense {ntW^ cm~^), embedded clouds (but see e.g. Aalto ct al. 1995 for 
a discussion on mid-IR pumping of HCN). The bulk of the molecular mass is 
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Figure 3. Cartoon models of a high pressure and low pressure ISM 

in the dense gas, but the high value of TZi-o can primarily be attributed to 
the large linewidths of the diffuse molecular gas (which can remain molecular 
because of the extreme pressure). Galaxies which may be characterised by this 
ISM scenario include Arp 220, NGC 6240 and Mrk 231. 

It may seem surprising that the ^^CO emitting gas is somewhat denser, 
10^ ^ 5 X 10^ cm~^, in the lower pressure scenario (and with a smaller volume 
filling factor). This may be because (in contrast to the "raisin roll" ISM) the 
intercloud medium is likely atomic (or ionized) since the destructive forces of 
the newborn stars are not balanced by a high ambient pressure. Thus, the ^^CO 
emitting gas is more "cloudy" in the lower pressure model, but the relatively 
faint HCN emission indicates that the mass fraction of high density gas, n > 10^ 
cm"^, is considerably lower. The ^^CO emission may, however, be filamentary 
and not always in bound clouds. The large value of TZi-o is likely caused by 
elevated gas temperatures (Tk^SO K) since the gas is dense enough for ^^CO 
to be thermalized (Tex ~ Tk). The clouds may be PDRs (photon dominated 
regions) if the impact of the starburst is strong enough. Galaxies which may be 
characterised by this ISM scenario include the Medusa (NGC 4194), UGC 2866 
and NGC 1614. 



3. HCN — tracing star formation or pressure? 

Because of its high dipole moment, HCN 1-0 emission is often used as a tracer of 
high density (n > 10^ cm~^) gas and the line intensity ratio as a measure 

of the mass fraction of dense gas. Bright HCN (low ratios) in ULIRGs 

has been used to argue for a starburst origin of their luminosity (e.g. Solomon, 
Downes and Radford 1992) because of a global correlation between Lhcn and 
LpiR. It is therefore interesting that the SFE of the Medusa rivals that of the 
ULIRG Arp 220 — even though the HCN emission is faint towards NGC 4194. 
How can there be an efficient transformation of gas into stars when there is only 
little dense gas present? The answer may lie in the time the molecular gas spends 
in a dense phase which is related to the average gas pressure and dynamics. The 
Medusa starburst takes place in an environment of reduced shear because it 
occurs inside the region of solid body rotation. Thus, gravitational instability 
may dominate over tidal shear with a resulting increase in the SFE. The rate 
of star formation (SFR) is indeed very high, 40 Mq yr~^, which is close to 
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the maximum SFR per kpc^ found by Lehnert & Heckman (1996). Since the 
hydrostatic pressure is much lower than in a ULIRG, gas is not maintained at 
high densities, but passes through a dense phase quickly on its way to becoming 
stars. The time spent in a dense phase is short and will not be an observed 
signature of the ISM. The rate of star formation has been found to be dependent 
on large scale dynamics (e.g. Kenney, Carlstrom and Young 1993). Perhaps, for 
a threshold level of the gas surface density, the dynamical environment of the 
clouds is also important in regulating the efficiency of star formation, even more 
so than the relative amount of high density gas. 

For the compact burst in the deep potential well of a massive major merger, 
the gas surface density (and matter density) is high and thus also the ambient 
pressure leading to large average gas densities. Still, the dynamical and radiative 
environment may be unfavourable to star formation. Within 200 pc of the deep 
potential well of a differentially rotating galaxy with rotational speed 250 km s~^ , 
the clouds must have average densities of n > 10^ cm~^ just to be stable against 
tidal shearing. Lower density gas will only exist as diffuse unbound clouds. 
In addition, the feedback mechanisms of the starburst itself may help regulate 
the SFE, and these mechanisms may be effective in the densely packed, high 
gas surface density central regions of ULIRGs. Thus, in some circumstances, 
the gas will remain at substantial average density without forming stars. In 
this context, the i^HCN - -^fir correlation is, at least partially, caused by high 
pressure gas being more centrally concentrated than low pressure gas, and that 
FIR emission in galaxies tend to emerge from the inner regions — at least the 
60 fxm emission. The FIR emission may indeed originate in starburst activity, 
but also from an AGN or from an evolved starburst where densely packed stars 
heat a fragmented ISM. 

A possible caveat in the notion that the SFE is similar in the Medusa and 
Arp 220 is the possibility that the -^^CO luminosity is tracing molecular mass 
differently in a compact and extended starburst. Furthermore, the ULIRGs 
deviate from the limiting SFR per kpc^ found by Lehnert and Heckman. Their 
maximum star formation rates are in the range 100-300 M© yr"^ instead of 
20-40 — and from a smaller area than less powerful galaxies. Dust opacity and 
hidden AGNs may explain some of the discrepancy in the limiting SFRs, but we 
cannot exclude that the underlying star formation mechanisms are different. 
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